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Abstract 

It is shown that analytic properties of standard QCD perturbation 
theory contradict known spectral properties, and contain in particu- 
lar IR generated ghost poles and cuts. As an outcome the rigorous 
background perturbation theory is developed and its analytic prop- 
erties are shown to be in agreement with general requirements. In 
a limiting case of large Ac, when QCD amplitudes contain only pole 
singularities, the strong coupling constant as{Q) is shown to be also 
meromorfic function of external momenta. 

Some simple models and examples are given when nonperturba- 
tive (5 function and as{Q) can be written explicitly. General form of 
amplitudes at large Ac is given in the framework of background per- 
turbation theory and its correspondence with standard perturbation 
theory at large momenta is demonstrated in the example of e^e" an- 
nihilation. For time-like momenta the background coupling constant 
differs drastically from the standard one but the background series av- 
eraged over energy intervals has the same (AF) asymptotics at large 
momenta as in the standard perturbation theory. 



1 Introduction 

The standard perturbation theory (SPT) in QCD is well developed both 
on the theoretical and on phenomenological levels |31;0 ? and constitutes 
the major and best understood part of QCD. Successful applications of per- 
turbative expansions for processes with large momenta are numerous and 
impressive. 
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Howerever there are several basic difficulties in SPT (it is assumed that 
RG improvements such as partial summation of large logarithms are auto- 
matically included in SPT), namely 

1) Analytic properties of SPT amplitudes do not correspond to the ex- 
pected spectral behaviour. In particular, there appear ghost poles and cuts 
even in the Euclidean region of momenta, where one expects amplitudes to 
be holomorphic. 

These ghost singularities are due to analytic properties of RG improved 
as{Q) and are connected with IR divergence of as{Q)- Also in the Minkowskian 
region of momenta as{Q) have logarithmic cuts which the physical amplitude 
should not possess, and they are therefore artefact of SPT. 

In short, amplitudes computed in SPT have analytic behaviour which has 
nothing to do with physical thresholds and cuts due to creation of hadrons. 
Therefore one may only speak about some duality relations between SPT 
amplitudes and physical amplitudes integrated over sufficiently large energy 
interval. 

2) Another basic defect of SPT is the lack of convergence of standard 
perturbative series. There are arguments that the latter is an asymptotic 
series 0, and in all cases it is not clear where the series should be cut offi 
(Some hints that the three-loop contribution deteriorates the physical results 
in the 1-1.5 GeV region are contained in and [Q). Moreover, the Landau 
ghost poles give rise to the appearence of the so-called IR renormalons 0, 
which make the sum of perturbative series undefined even in the Borel sense. 

The attempts to connect the IR renormalons with the nonperturbative 
contributions may have only qualitative character at best. Therefore strictly 
speaking the notion of the sum of SPT series has no definite meaning, and 
one may only hope that at large Euclidean momenta the first few terms of RG 
improved perturbative series describe the asymptotics of physical amplitudes 
with reasonable accuracy. 

To improve the situation several approaches have been developed. First 
of all, it was understood long ago, that in addition to SPT also nonperturba- 
tive contributions should be taken into account. Technically the latter have 
been introduced as the local terms (condensates) in the OPE, and on this 
foundation the QCD sum rules have been introduced ||^. In this way one can 
approach the low-energy region around 1-2 GeV and a lot of useful physical 
information have been obtained in this method for the last 20 years 0. 

However in this method one does not solve both problems of SPT de- 
scribed above, and instead postulates that in addition to a few ffist SPT 



2 



terms one can add few power terms to imitate behaviour of physical amph- 
tude not only in asymptotic region but also in the region of few GeV. In the 
lower energy region the OPE series has at least double divergence: due to 
explosion of as near the Landau ghost pole and due to explosion of power 
terms 

It is important to note that the original OPE and QCD sum rules Q have 
been properly defined only in the Euclidean region, and the transition to the 
time-like region is assumed to be done aposteriori, after all calculations are 
done in the Eulidean region. 

Definition of as{Q'^) in the Minkowskian region of Q^, < 0, is a problem 
by itself in SPT, since formal analytic continuation of two - and more-loop 
expressions for ^^(Q^) yield complex expressions violating exphcit unitarity 
conditions in nonasymptotic region. 

Moreover as was stressed in 0, the invariant coupling g{Q) can be defined 
only in the space-like domain, and "inside the RG formalism there is no 
simple means for defining g{Q) in the time-like region". 

New formalism for definition of as(Q^) both in space-like and time-like 
region was developed in |Q and a new method was suggested, called analytic 
perturbation theory (see ^ for a review and further references), where q;s(Q^) 
is forced to be analytic for > and a special procedure is envisaged to 
continue ag analytically into the time-like region. The method allows to get 
rid of the Landau ghost pole and cuts and seriously improves convergence of 
the perturbative series. 

It is clear however that this is not unique way of analytic definition in 
the whole complex plane of z. 

In what follows we shall choose a completely different strategy. To sim- 
plify matter we shall consider below the limiting case of large N^. In this 
case one can be sure that all physical amplitudes contain only poles as func- 
tions of external momenta |jlO|], and we shall require that perturbative expan- 
sion would reproduce the meromorphic properties of physical amplitudes, i.e. 
as(Q^) have only singularities in the time-like region, and those are poles. 

To make this goal, one needs to use Background Perturbation Theory 
(BPT) instead of the standard one, and we shall derive rigorous formalism 
based on BPT developed in the seventies and eighties |]TT| and generalized in 



12|-|14| to include nonclassical background and averaging over background 



fields. 

In this way one obtains a systematic formalism which allows to express 
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all terms of BPT through the irreducible correlators of background fields 
(integrals thereof) and renormalized coupling constant as, which we shall 
denote asis) to distinguish it from as{s) in SPT. 

It was argued in that asis) and physical amplitudes satisfies the same 
RG equations, in particular Ovsyannikov-Callan-Simanzyk (OCS) equations, 
and the important difference in a^lQ^) from as(Q^) lies in the character of 
its dependence on Q^. It was demonstrated in |jl2|-|I^ that asiQ^) has the 
property of freezing, or saturation at small Q^, i.e. it tends to a finite limit 
ctB(O) when — >■ and has no singularities in the whole Euclidean region 
Q2>o. 

This behaviour of a_B((5^) was tested repeatedly in e~^e~ annihilation 
in fine structure of charmonium and bottomonioum levels \TE 



and recently 0] in comparison with accurate lattice data on small-distance 
behaviour of aiiR)- In all cases the same form of solution for as was used 
without free parameters which produced results in a good agreement with 
experimental and lattice data. In this way the phenomenon of freezing (sat- 
uration) was demonstrated both theoretically as a result of confinement in 
background fields [0, and phenomenologically in comparison with lattice 
and experimental data. 

With all that some important theoretical questions were not answered, 
first of all what happens with asis) in the Minkowskian (time-like) region 
and what are analytic properties of asis) in the whole s-plane. 

Secondly, what is connection of asis) and other nonperturbative defini- 
tions of a{s), e.g. lattice definitions of a — aiis). 

In the present paper it is our purpose to study the problem of BPT in 
the large Nc limit, trying to elucidate several aspects. 

Firstly, we formulate the foundations of BPT and rules for calculation of 
perturbative series, starting from the purely nonperturbative term. 

Secondly, we derive the RG equations taking into account that nonper- 
turbative background is in general not classical solution and subject to the 
vacuum averaging, which due to the 'tHooft identity can be done indepen- 
dently of the perturbative field averaging. 

Thirdly, we find the most general solution of RG equations and in particu- 
lar the nonperturbative (3 - function which has known lowest term expansion, 
at the same time ^^(Q^) is represented as a sum of pole terms. 

To understand which kind of singularities asiQ"^) niay have compati- 
ble with analytic properties of physical amplitudes we consider a simplified 
model and demonstrate that any finite order perturbative expansion has ad- 
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ditional singularities which are ehminated when a partial summation of the 
perturbative series is done. 

As a next step we formulate the generic perturbative expansion for a phys- 
ical amplitude, choosing an example of e^e~ annihilation into hadrons and 
write it as a sum over poles in the time-like region with calculable coefficients. 

At this point one may wonder how this meromorphic expansion is con- 
nected to SPT expansion, where the perturbative series is in powers of 
as{Q'^), which contain logarithmic functions of and not meromorphic. 
The answer to this connection was given in the analysis of the lowest term of 
e'^e" annihilation - the hadronic part of the photon self-energy. The latter is 
a sum over meson poles at large Nf. and has a proper logarithmic behaviour 
logQ^ at large Q^. We demonstrate that a similar correspondence takes place 
also for higher terms of BPT, and formulate conditions on the coefficients 
of meromorphic expansion which are necessary to reproduce the known SPT 
expansion at large Q^. 

Discussion of practical applications and comparison to other approaches 
concludes the paper. 

The exact structure of paper is as follows. In section 2 general rules of 
BPT are given, and Green's functions for valence quarks, gluons and hadrons 
are written explicitly. 

In section 3 the RG equations are written down in the case of most general 
background and the point of RG scheme is discussed and a particular solution 
of RG equations for f3- function and ^^(Q^) is presented. 

In section 4 a simple model is discussed where some BPT subseries can 
be summed up exphcitly. 

In section 5 nonasymptotic terms in the spectral sums are calculated and 
in section 6 behaviour of perturbtive series in the region of time-like is 
studied in detail. 

Discussion and outlook is contained in the concluding section. 

2 Background Perturbation Theory for arbi- 
trary background 

Gluon field plays in QCD two different roles: 

i) gluons are propagating, and at small distances this process can be 
described perturbatively, leading in particular to color Coulomb interaction 
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between quarks (antiquarks); 

ii) gluons form a kind of condensate, which serves as a background for the 
propagating perturbative gluons and quarks. This background is Euchdean 
and ensures phenomena of confinement and CSB. 

Correspondingly we shall separate the total gluonic field into pertur- 
bative part and nonperturbative (NP) background B^: 

= 5^ + (1) 

There are many questions about this separation, which may be answered 
now only partially, e.g. what exactly is the criterion of separation. Possi- 
ble answer is that perturbative fields get their dimension from distance 
(momentum), and therefore all correlators of fields (in absence of B^) are 
singular and made of inverse powers of (x — y) and logarithms, where enters 
the only dimensional parameter of perturbative QCD - Aqcd- Therefore 
evidently any dimensionful constant, like hadronic masses or string tension 
cannot be obtained as a perturbation series. In contrast to that, NP fields 5^ 
have mass dimension due to the violation of scale invariance which is intrin- 
sically present in the gluodynamics Lagrangian. The origin of separation (|I|) 



is clearly seen in the solutions of nonlinear equations for field correlators IjT^] : 
a perturbative solution of those leads to singular power-like field correlator, 
whereas at large distances there is a selfconsistent solution of the equations, 
decaying exponentially with distance with arbitrary mass scale, since equa- 
tions in are scale-invariant. Full solution including intermediate dis- 



tances produces mixed perturbative-nonperturbative terms, containing both 
inverse powers of distance and exponentials. For these terms criterion of 
separation fails. 

One can avoid formally the question of separation principle (and of dou- 
ble counting) using t'Hooft identity [T^, which allows to integrate in (1) 



independently over 5^ and a^: 

Z = ^J DB^7]{B)DiPD^Da^e''^''^ (2) 

Here the weight ri{B) is arbitrary and may be taken as constant. 

To define the perturbation theory series in (^a^ one starts from (0) and 
rewrites the Lagrangian as follows: 

Ltot = Lgf + Lgh + L{B + a) = 
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Lq + Li + L2 + Lint + Lgf + Lgh (3) 
where Li have the form: 

= ^at[DTDl%,-D^;Df-grX^]a', ' (4) 

Lo = —{FHAB)? ; L, = alDl\B)F;, 
Unt = -l{D,{B)a, - D,{B)a,rg r'^a^al - -/ r'^'o'^airf a^at (5) 

Lghost = -e+{D^{B)D^{B + a))abeb. (6) 
It is convenient to prescribe to B^, the following gauge transformations 

a^, U^a^U (7) 

B^ ^ U{B^ + '-d^)U (8) 

and to impose on the background gauge condition 

^ {D,a,r = d,al + gr^^Bla^ = (9) 

In this case the ghost field have to be introduced as in and the gauge- 
fixing term is Lgf = — ^(C)^ One can write the resulting partition function 
as 

Z = ^lDBrjiB)ZiJ,B) (10) 

where 

Z{J,B) = J DipDipDa^DeDe+e^""^!-^^''^'^'' (11) 

We now can identify the propagator of from the quadratic terms in 
Lagrangian ^2(0) — ^(G")^ 

Gti = [^A V - D,D, + igP,, + -^b^D,]-J (12) 
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It will be convenient sometimes to choose ^ = 1 and end up with the well- 
known form of propagator in - what one would call - the background Feyn- 
man gauge 

G''X = {Dl-5,.~2igF,,)-' (13) 
Integration over ghost and gluon degrees of freedom in ([Til) yields 
Z(J, B) = const{detW{B));^J^[det{-D^{B)D^{B + a)]^^x x (14) 



^{^+j:^ia = jj)}exp{-^JW-'j) 



J.=D^{B)F^,{B) 



where W = G-\ and Gis defined in 

Let us mention the important property of the background Lagrangian 
(|^): under gauge tranformations the fields a^^B^ transform as in (0-^ and 
all terms of (|^), including the gauge fixing one |(G")^ are gauge invariant. 
That was actually one of the aims put forward by G. 'tHooft in |Tl|, and it 



has important consequences too: (i) any amplitude in perturbative expansion 
in ga^ of ([TT|) and (p!4|) corresponding to a generalized Feynman diagram, is 
separately gauge invariant (for colorless initial and final states of course), 
(ii) Due to gauge invariance of all terms, the renormalization is specifically 



simple in the background field formalism |rT[, since the counterterms enter 
only in gauge-invariant combination, e.g. and hence Z-factors Zg and 
Za are connected. We shall exploit this fact in Sections 6 and 7. 

Let us turn now to the term Li in (|^) Li = 2tr{ayD^F^y), which is 
usually missing in the standard background field formalism 



tacitly assumes there that the background B^ is a classical solution, 

D^{B)F^,{B) = 0. (15) 

Here we do not impose the condition (|T5|) and consider any background, 
classical or purely quantum fluctuations. Let us estimate the influence of the 
vertex L^. In general it leads to the shift of the current in the expression 
for the perturbative series ([M])- Physically it means that at each point the 
background B^ can generate a perturbative gluon via the vertex {a^DyFi,^), 
and this vertex is proportional to the degree of " nonclassicality" of B^. For 
the quasiclassical vacuum, like the instanton model, the average value 
< {D^F^^Y > over instanton ensemble is less then Q{p^ / R^) and is small ( 
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at most of the order of few percents) while average of < D^F^^ > vanishes in 
the symmetric vacuum. All this is true provided the instanton gas stabilizes 
at small density. 

Let us estimate the effect of Li in the general quantum case. To this end 
we calculate as in [12] the contribution of Li to the gluon propagator. 

If one denotes by the integral Da^ with the weight Lfot as in (|ll]), 
we obtain 

< af,,{xi)a^^{x2) >a= G^,^^{xi,X2) + A^^^^ 



J d'^yid'^y2Gf,^y^{xi,yi)DpFp^^{yi)DxFxy^{y2)Gy^f,^{y2,X2) (16) 



The gluon Green's function G^^ is given in ([T3|) and depends on the back- 
ground field Bp, as well as Dp and Fpx. To get a simple estimate of A we 
replace Gp^ by the free Green's function G^p} and take into account that 

< DpFp^^{yi)DxFxu,{y2) >b^ — < Fp^^{yi)Fxu,{y2) > (17) 

oyip oy2\ 

and for the latter we use the representation |TB[ in terms of two independent 
Lorentz structures, Dijji — 1/2) and -Di(i/i — 2/2)- The contribution of D (that 
of D\ is of similar character) in the momentum space is 



AS,.(fc) - '-^^ll^l^^J^D{k) (18) 

where 

D{k) = J d'ye''yD{y) (19) 

Insertion in (|19|)of the exponential fall-off for D{y) found in lattice calcula- 
tions [|19] yields 

With /i ^ IGeV, D{k) ^ 0.12. Thus A'^{k) is a soft correction fast 
decreasing with k as k~^ to the perturbative gluon propagator. 
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Renormalization group improvement of per- 
turbative series and analytic properties of 



The most important property of the Background method discussed in the 
previous chapter, is the gauge invariance of the total Lagrangian including 
the gauge-fixing term Lgf. 

As a consequence of this the counterterms in the renormalization pro- 
cedure have to be gauge invariant too and this fact establishes connection 
between Z-factors of charge and background field |Tl| , namely, if 



go = Zgg, Bf = Z]i'B^ (21) 

then 

ZgZ]i^ = 1. (22) 

This property is very important for the BPT and the Field Correlator 
Method in general. Indeed the combinations entering in the latter have the 
form 



A„ = tr(^F(S(l)<|.(l,2)^F(f)^(2)<|.(2,3)...^F^^^^ (23) 

where Fj^^^^ is made of the background field B^ only, and parallel transporters 
$(i, fc) making Eq.(p3D gauge invariant depend also on gB^. Hence all ex- 
pression in (^Bp contains only combinations gB^, which according to (P^) are 
renorminvar iant . 

Therefore in all RG relations, in particular in Ovsyannikov-Callan-Simanzyk 
(OCS) equations the background fields (averaged as in A„ or nonaveraged) 
enter on the same ground as the external momenta, and hence all relations 
can be kept intact. 



E.g. OCS equations for some physical amplitude F now look like |14 



+ (3— + m(7^ - 1)- 7r I F(Api, ...Apjv, A^'^A^, g,v) = Q 



\ d log A dg"^ dm 

(24) 

where we have defined as usual 
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and similarly for 7^, with boundary conditions 

a,(A = 1) = a,{v). (26) 



In writing the r.h.s. of (^q) we tacitly assume that depends also on all 
external parameters {pi}, {A^} in addition to the dimensional regularization 
parameter v (or its equivalent in other schemes /x). The reason why exter- 
nal parameters are not usually written lies in the fact that RG equations 



(like (p4D, (p5|)) define only dependence on one scale parameter, e.g. A and 
through it on v (/x), while dependence on external parameters is not fixed, 
and they may in principle enter in arbitrary dimensionless combinations, 
unless additional information is obtained from perturbation expansion (BPT 
or SPT). 

Several words should be said about the choice of renormalization scheme 
for BPT. 

It is convenient to use the MS scheme, since also in the presence of 
background the modified Feynman amplitudes of BPT are diverging at small 
distances (large momenta) and, when these momenta are much larger than 
the average magnitude of the background fields, one can neglect the latter, 
e.g. using OPE one expects corrections of the form — which do not alter 
Z factors. In particular the difference F — Tpert = AF, where Tpert is a usual 
sum of Feynman amplitudes without background, is UV convergent, so that 
all UV divergencies of F are the same as in Tpert- 

Now we turn to the possible solutions of Eq. (pSf). The SPT expansion 
for P{as) looks hke 

^ ("^^ = -2;^"^ - 4^"^ - 64^"^ - (4^ - - ^'^^ 
where the first coefficients areQ 



5f/(3); Uf = 0; bo = 11; h = 51, 62(M5) = 2857, h{MS) = 58486. (28) 

sum- rif = 0- bo = yiV,; 6i = ^Nl (29) 
Following (^) it is reasonable to look for (3{as) in the form 



^Note that definition of our coefficients bi differ from that in ^ and coincide with [30] 
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where (f{x) is an unknown function and <f'{x) is the derivative in the argu- 
ment X. Solving (|30|) for one obtains 

4:71 

= MWlWfRXj] (31) 

where we have replaced X ^, and is an arbitrary function of external 
parameters of dimension (mass)~^. 

/^From p7| ) one concludes that the only condition on (p is the expansion 



in powers of 



One can rewrite (|3lD in the form 



X(Q') = ^ - ^V^(^) (33) 

where z = 1/as and x(Q^) = 1^ In/x^C^, and is the representative of the 
external parameters. 

As a simple example one can consider the self-energy photon function 
n(Q^) entering in the e~^e~ annihilation process, in which case is the only 
external momentum, while background fields create masses of the bound 
state spectrum. In this case one can take the form appearing in the lowest 



approximation of BPT |jTj], leading to the replacement 




(34) 



Here ■ip{x) = is the Euler function, having only simple poles at 

X = —n, n = 0, 1, 2, .... 

This example clearly shows that x(Q^) is a meromorphic function of Q^, 
and we must require that as be also a meromorphic function of Q^, analytic 
in the Euclidean region ReQ^ > 0, since otherwise the BPT would violate 
the analytic properties of amplitudes in the limit Nc oo. 

In this way we are lead to the conclusion that <f{z) should be a mero- 
morphic function of 2; = since in this case the r.h.s. of Eq.(P^ is a mero- 
morphic function of (meromorphic function of a meromorphic function of 
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some argument is again meromorphic) . Thus one can equalize functions on 
both sides of Eq.(0) since both have the same analytic properties. 

We do not specify at this point the character of possible singularities; the 
resulting physical amplitude should have only simple poles at bound state 
energies, while any finite series of BPT may have poles of higher order and 
additional poles, as will be seen in examples below. 

To illustrate possible forms of ^p{z) let us consider one example of mero- 
morphic function, namely 



(fi (—^ = ip (— + aV a = const. 



Perturbative expansion of ip has the form 

^P'{z) 

or applying it to (^) one has 



1 1 1 

z 2z'^ 



^' ( - + A 



a. 



z —y oo 



a" 



l + a,A 2(l + a,A)2 6(1 + a,A)^ 



+ 



(35) 



(36) 



(37) 



Therefore one obtains 
1 



i-^'(^ + A) 

■bo ^ ^as > 



1 + ^r—;-(^s + 



27r6n 



-A + - + 

2 27r6o 



+ 



6i o 



A^ - A + - + 
6 



27r6n 



;i-2A) 



27r6n 



(38) 



^From (|38|) one can find the "theoretical" value of 62 and 63, namely 



= 647r26i[A2- A + ^ + 27 



3173.4 - 167r6iA 



(39) 



A +7^ 
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27r6n 



Now since '?/'(z) has poles at z = 0, —1, —2, + A^ has poles at the 
values of as 

"^ = -1'-TTa'-2Ta'- 
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Correspondingly ip'i^) has poles of second order, and as a consequence of 
(|30| ) P{as) has zeros at the values of given in (1^). These zeros condense 
to tts = from the negative side. 

Poles of P{as) are defined by the zeros of the denominator in (pOD, i.e. 
given by 

1 = Ai^' f- + aV (41) 



27r6o \«„ 

For A > and > there is one root of (|TD at as = a;^^\ if A < Aq 
where Aq is defined by 

V''(Ao) = — ^ = 1.355; Aq = 1.145. (42) 
bi 

For A < Aq one obtains one pole of P{as) given by (^TD, while for A > Aq 
there are no poles of P{as) for > 0, and P{as) is monotonically decreasing 
function. 

For as < one has infinite number of poles of f3{as) which lie between 
zeros of /3{as), described by Eq. (^Oj). 

We now come to the analytic properties of as{Q^), having in mind parametriza- 
tion (|3^ ) for Il{Q^) in e+e~ annihilation. 

Poles of as are to be found from the equation 

and we require that for ReQ^ > there should be no poles of asiQ"^), which 
gives condition 

Now ij{z) > for ;z > 1.46, and for = = lGeV^;A = 0.37^6^ 
(the accurate calculations in [|, 0, |1^ strongly prefer Mq = 1 GeV) one 
obtains for the sum of two first terms in (E3f) the value 



In = 1.99 - 0.577 = 1.41 (45) 

A? \im? ) ^ ' 

Finally one can see that for A > Aq one has > —0.5 and the l.h.s. 

of (|4^ is positive, and hence for all > also l.h.s. of Eq.(^3D is always 
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positive. In this way we have proved that <ys{Q^) has no poles for all > 
if A > Aq. 

The resulting form of OsiQ"^) valid in all plane is now 

«B = — ; 77W7172-, -, — ^ (46) 



lnf + *(2iS!Si) + tV'(^ + A 



where A > Aq = 1.145. For < the function as{Q ) in (|6|) has zero's 
and poles which not necessarily coincide with poles of ip ( — ^^^j-'^ j . 

But this is as it should be, since the amplitude Il{Q^) must have poles at 
the final positions 

-Q2 = (m("))2, n = 0,1,2,... (47) 
which take into account valence (and Coulomb) gluon effects. 

4 A simple model for constructing perturba- 
tion series 

In this section we shall consider a simple model which is to illustrate the 
analytic properties of separate terms of BPT and of its infinite sum. Our 



discussion in this section will be formally similar to the lectures pO[ where 
old-fashioned perturbation theory was used for SPT, except that we use 
confined states instead of free quarks and gluons. To this end we define a 
quantum-mechanical Hamiltonian Hq which produces bound states (equiva- 
lent of infinite number of mesonic states) and two second-quantized operators 
Vi(t) and V2{t) which create additional valence gluons, to be bound together 
with quark and antiquark in a sequence of hybrid states. We define the total 
Hamiltonian 

H = Ho + V,{t) + V2{t) (48) 

and 

Vi{t) = e'^'^'via^ + h.c. (49) 

V2{t) = e'^^^^b-^a+a + h.c. (50) 

where a, and 6, fo"*" are annihilation (creation operators). We do not specify 
the (i-space dependence of Vi{x),V2{x)\t- is the common time of the instant 
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form of dynamics. It is clear that Vi, V2 are prototypes of the terms ipdip and 
a^da in the QCD Lagrangian, here to simphfy matter we disregard in a^(x) 
or A/i(k) dependence on polarization fi and momentum k, but introduced 
b, a to distinguish gluon emission by quarks (in Vi) and by gluons (in V2). 

We assume that Hq has only discrete spectrum of mesons (confinement) , 
so for the Green's function we have 

z^Goit) = H,{t)Go; Go(t) = E^i^Vf ^e"^^^"^* (51) 
while for the Fourier transform one has 

GoiE) = / e-^^^Gom = z ^ ^t^' ^^2) 

"^0 n ti — tin 

For the full Green's function one can write 

i^G{t) = {H^ + V{t))G{t), V{t) = V^{t) + V2{t) (53) 
at 

with the solution 

G{q,,q2,t) gi|e-/o(^°+^(*'))'^*'|g, > . (54) 
The perturbative series looks like 

Git) -^Go + j Go{t - ti){-i)V{ti)Go{ti)dti+ 

+ 1 lGo{t-ti){-i)V{ti)dtiGo{ti-t2){-i)V{t2)Go{t2)dt2 + ... (55) 

For the energy- dependent Green's function G(E) the perturbative series 
can be written taking into account the vacuum state for a, b operators 

a\vac >— b\vac >— (56) 

which yields 

G{E) = Go{E) + Go{E){-iVi)Go{E - u;i){-iVi)Go{E) 

+G'o(^)(-^V\)G'o(^-c^i)(-^K2)G'o(S-cc;i-cc;2)(-^V^2)G'o(^-^i)(-^V^i)G'o(^)+ 
+Go{E){-iVi)Go{E - u;i){-iVi)Go{E){-zVi)Go{E - u;,){-iVi)Go{E)+ 

16 



+Go(£;)(-^V^i)Go(^-a;i)(-«ri)Go(£;-2a;i)(-«^i)Go(£;-cc.i)(-^V^i)Go(£;)+... 

(57) 

If one introduces the "gluon propagator" 

Ki{E) = {-iVi)G^{E - uj{){-iVi) (58) 

then the partial summation of diagrams with gluons emitted by quarks only 
yields 

G(^^)(E) = G,{E) + G,{E)K,{E)G,{E) + ... = Go{E)- (59) 

1 - AiGo 

Similarly for the "gluon-b" propagator 

K2{E) = {-iV2)Go{E -uoi- uJ2){-iV2) (60) 
and the "full gluon propagator" 

= K^{E) + {-iVi)Gq{E - uj^)[l + K^G^iE - cui)+ 

{K^GoiE - uj,)f + ...]{-iV,) = {-iVi)G,{E - u;,)- -A-iV,) 

(61) 

One should note that Vi are operators and denominators in (^) and (|6TD 
are formal operator expressions. Now to analyze the perturbative expan- 
sion (^) or its partial sums (^),(0) one should specify the unperturbed 
spectrum (0) and define matrix elements of Vi between unperturbed eigen- 
functions '^^n' ■ We shall assume for simplicity that that the spectrum is linear 
in n, 

£;(o) = Mo + mn, n = 0, 1, ... (62) 

and for V^^ = {ip'^'^^^ViLpf'^) one can assume that at large n, /c this matrix 
element factorizes, vj:^^ = c^'^^^c^*''. 
As a consequence one has 

where we have defined 

n,iE-u;,)^-Y. j, ''f'' , (64) 
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It is clear that Qi{E — ui) has poles at the "gluon excited" i.e. "hybrid 
state" masses E = Ek + cui, however the partial sum G^^'^\E) has shifted 
poles due to vanishing of denominator in (|63|) , in addition to unperturbed 
poles at En- 

A similar situation occurs for G^^^\E). If one defines 

n,{E-co,) = j: j, 'p' (66) 

U,{E -u:,-u,) = -Y. (67) 
then one can rewrite (^) as 

r(^2)/n _ i-'^Vl)\'Pn > rr I ^+(2) n2(^-^i -a;2)c(g^ ^ < V^m | ( 

E-E„-cu/ "^"^ " l-fi2(^-^i)n2(^-a;i-cu2)^E-E„ 

(68) 

This exercise was undertaken to illustrate that the "full gluon propagator" 
G^^^\E) has in addition to unperturbed "one-gluon poles" also 
shifted "two-gluon poles" due to vanishing of denominator in (0), but not 
two-gluon poles" themselves. 

To make our model more realistic and finally to incorporate in it QCD 
features, one should make several steps. Firstly, one must take into account 
negative energy states, corresponding to the background in time motion, i.e. 
as in the Feynman propagator we replace (cf.[pO|) 



1 1 1 



E-En E^-El s-iMi'Y 
Correspondingly, the hybrid propagators take the form 
1 1 1 



(69) 



E-Ui-En S-{M^nY E -UJ1-UJ2- En S-{M^n 



(2) A 2 



(70) 



where M'^\M'^^ are hybrid masses with one and two gluons respectively. 

Now one must take into account the structure of Vi? ^ operators in QCD. 
Both describe emission of one gluon from quark and gluon line respectively. 
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and both are local field operators, whereas wave functions refer to the nonlo- 
cal instantaneous bound state objects, namely, mesons (gg), 1-gluon hybrid 
{qqg) and 2-gluon hybrid (qqgg). We shall denote the full set of the corre- 
sponding eigenfunctions as ip!^\x), ip!^J^^{x, y) and ^n)^^u2(^^ Note that 
the string pieces between quarks and gluons formed by background field are 
straight fines, given by vectors x, y in (p!^\ and by x, y, z in (p!^\ and they 
contribute to the rotational and vibrational eigenenergies. We also note that 
the string vibration in this background picture is described by the full set of 
one-, two- and many-gluon hybrid states. 

Here we have introduced in hybrid eigenfunctions additional lower in- 
dices i/i and i'i,i'2 to describe the quantum numbers of additional bound 
gluons, note that i^i are multidimensional vectors with integer components. 
Correspondingly hybrid masses depend on those quantum numbers; 

M(o) = M(°)K,L;J^^) (71) 

Mi'^^M^'\nr,L-X,L';J^^) (72) 

Mi'^^M^'\nr,LX,L',n';,L";J^^). (73) 

Here UrL refer to the radial quantum number and angular momentum of 
quark with respect to antiquark, each additional gluon has its own n^, L with 
respect to neighboring quark or gluon. We have not taken into account yet 
spin degrees of freedom of q, q and gluons, which create fine and hyperfine 
structure of levels, and we shall neglect those, since we are interested in the 
dominant contributions to the masses at large excitations, and our final goal 
is the result of summation over all excitations in G{E) and establishing of 
some relation between SPT and BPT, which requires analysis of high excited 
states in the sums. 

Now the local character of the operators Vi,V2 mentioned above leads 
to the following properties of matrix elements (since Vi acts only at one of 
the end points of vectors xy, it does not influence quantum wave function, 
unless Vi contains derivative in coordinate as in V2, and in this case it changes 
quantum numbers by a discrete number): 

< nr,L\V^^^\nr,L;n'^,L' >~ SnrUr^LL 

< Ur, L;nl, L'\V^^')\nr,L,n'^,L';n';, L" >~ ^n^n^^LL^^^i'L'- (74) 

Moreover we shall assume that {M^^y depends on all quantum numbers 
linearly, 

(M«)2 = m2(a^n« + 65]L«) (75) 
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Let us now apply these rules to the computation of the "full gluon propaga- 
tor" and try to find a correspondence with the standard perturbative 
calculations. 

To this end we assume that the sum entering in Ki{E) is converging, i.e. 
/ (T^w ^ ^<n\-iVi\n,v><n,v\-iVi\m> 1 

<u\ME)\^>=Y. grrWiKJ) 

where we denote n(i/) = n^, L{n'^, V) and replace s —Q^. 
The same type of expression for K2{E) yields 

ir^/rnM- ~ <n,v\-iV2\n,v,v' ><n,i>,v'\-iV2\n,i) > 
< n, .IK.mn, . >= E, + ' 

(77) 

Due to the properties (|75D the sum in (|77D is actually a double sum in, say, 
L'j,,?T,^, while depends on them linearly and we assume at this point that 
coefficients in ( ffTl ) do not depend on v' . Therefore one must regularize this 
expression, making one subtraction at s = 0, and renormalize the retaining 
logarithmic divergence. This is in full correspondence with the standard 
perturbative calculation of one gluon loop correction to gluon propagator. 
As a result K^iE) aquires the form 



< n, v\k2{E)\n, V >~ In { ^'^ ^'^^^ ^^' ) 6^J^, (78) 



where 

SM' = Ml^{n,i.,u')-M'in,i.) (79) 

and minimization is taken with respect to u'. 

Now one can consider G2{E) extracting coupling constant Qq from Vi, V2 
one has a series (at large s) 

Expression ( ppj ) is similar to that of the running coupling constant, but with 
the mass M^, corresponding to the freezing coupling constant derived in 
||r2| and studied later in [|l^, |T^. We note that logarithms in (^) appear 
as a result of summation over the string like spectrum (^) with the coef- 
ficients asymptotically constant, which is the property of matrix elements 
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of eigenfunctions of linear interaction [Q. Therefore this property and the 
correspondence found above can be only an asymptotics, while at finite s 
there are corrections, which should be taken into account. In the following 
chapter we look into this problem more carefully. 



5 Nonasymptotic terms in the spectral sums 

We start with the definition of the photon vacuum polarization function 
n(g^) and after a short review of asymptotic correspondence between the 
spectral sum and perturbative expansion, we shall look more carefully into 
the correction terms and comparison between OPE and spectral sum. 

We define the photon vacuum polarization function Il{q^) as a correlator 
of electromagnetic currents for the process e~^e~ hadrons in usual way 

- 2 J d^xe^i^ < 0|T(j^(x)j,(0)|0 >= {q,q, - g ,uq^)Yi{q^) , (81) 
where the imaginary part of 11 is related to the total hadronic ratio R as 

= -'f; - ""f "f = 12.e/mn(4. aM) (82) 

There are two standard approaches to calculate n(g^); first, when purely 
perturbative expansion is used, which is now known to the order O(a^) [21]. 
The second one is the OPE approach [6] , which includes the NP contributions 
in the form of local condensates. For two light quarks of equal masses (m„ = 
nid = m) it yields for 11 (Q^) 

;|n(g) = -^(l + -)/n- + ^+ + ^^^^^ +... (83) 

To make explicit the renormalization of a^, one can write the perturbative 
expansion of the function Il{Q^) (23) as 

n(g2) = n(o)(g2) + aBii^'\Q') + alu^^\Q') + ... (84) 

We now again use the large Nc approximation, in which case 

n(o) 

contains 

only simple poles in [27,28]: 
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and the mass M„ is an eigenvalue of the Hamiltonian H^^\ which contains 
only quarks and background field 5^, 

= M„^„, (86) 

In what follows we are mostly interested in the long-distance effective 
Hamiltonian which obtains from Ggg for large distances, r ^ Tg, where Tg is 
the gluonic correlation length of the vacuum, Tg a; 0.2/m, [19]: 

Gggix,0) = I DBr]{B)tri^^Gg{x,Oh^GgiO,x)) =< x|e-^*"'l"l |0 > (87) 

At these distances one can neglect in Gq the quark spin insertions a^^F^^ 
and use the area law: 

< Wc >^ exp{-aS„,in) (88) 

where Smin is the minimal area inside the loop C. 

Then the Hamiltonian in (|87| ) is readily obtained by the method of [22]. 
In the cm. system for the orbital momentum / = it has the familiar form: 

H^'^^ = 2^f + m? + err + const, (89) 

where constant appears due to the self-energy parts of quarks [23], while 
for / = 2 a small correction from the rotating string appears [22], which we 
neglect in the first approximation. 

Now we can use the results of the quasiclassical analysis of H^'^'> [24] where 
the values of M„, C„ have been already found. 

These results can be represented as follows (n = nr+l/2, = 0, 1, 2, / = 
0,2), 

Ml = 2na{2nr + I) + M^, (90) 

where Mq is a weak function of quantum numbers n^, / separately, comprising 
the constant term of (^); in what follows we shall put it equal to the p-meson 
mass, Mq ^ (see [23] for more discussion). 
For Cn one obtains quasiclassically [24] 

Cn{l = 0) = '^elN,m\ rr? = ina, 

Cn{l = 2) = ^elN,m\ (91) 
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Using asymptotic expressions (PI1|) , (PT|) for M„, C„ and starting with n = uq, 
one can write 

+12^2A.M2 + g2 12,^2'^^ m2 ^ ^^^^ 

+ divergent constant, 
where we have used the equahty 

nt;JoK^ + ^?^ ^ ^ ^ ^ 

+ divergent constant 

and = 

In (p2D we have separated the first no terms to treat them nonquasiclassi- 
cally, while keeping for other states with n > ng the quasiclassical expressions 
(pop, (pip. In what follows, however, we shall put no = 1 for simplicity. It was 
shown in [14] that even in this case our results reproduce e+e~ experimental 
data with good accuracy. 

Consider now the asymptotics of U^^^Q"^) at large Q^. Using the asymp- 
totics of ip{z): 



ij{z)z^oo = lnz---Y.^ ^ (94) 



where B„ are Bernoulli numbers, from (021) one obtains 



One can easily see that this term coincides at ^ Mq with the first 
term in OPE (22) - the logarithmic one. Taking the imaginary part of (39) 
at — > —s one finds 

R{q^) = 127r/mn(°)(-s) = iV^eg, (96) 

i.e. it means that we have got from U^^^ the same result as for the free quarks. 
This fact is the explicit manifestation of the quark-hadron duality. 
We can now identify U^^^Q"^) as 
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and compare this expression with the OPE series 

To this end one can use the asymptotic expansion of ip{x) 



'ipiz),^oo = In^; - — - ^ 



1 ^^2k 
k=l 



2z 2kz^^ 



1111 
In^ - + - + ... 98) 



and the expansion of 



oo 



*'W = -C-E(ri^-^). C = 0.577. (99) 



k=0 



Two features are immediately seen when comparing expansion (^) with 
the OPE (|83D (see also discussion in O, M) i) for m = the term 0(1/(5^) 



is absent, since one cannot construct a local gauge-invariant operator of di- 
mension mass squared. One should stress, that this type of operator may 
appear from the interference of perturbative and nonperturbative contribu- 



tions and is welcome on phenomenological grounds |2^ . However in the 
framework of BPT accepted in this paper the term U^^^Q"^) is solely due to 
NP contributions, and cannot have OPE terms 0{1/Q^). 

ii) The generic terms in Yl^^^Q"^) from expansion (R^) have the magnitude 

2ri 

^2Tr, while the corresponding OPE terms are much smaller, e.g. one should 
compare -f ^ (OPE) with ^ (spectral 

decomposition (0)). 

In the remaining part of this section we shall discuss possible solutions of 
problems i) and ii). 

Concerning the point i), one should say that (^7^ is an approximate ex- 
pression for H^^^ valid for large Q^, i.e. reproducing correctly the logarithmic 
term in (|95|) . However, if one interested in the next asymptotic terms, one 
should take into account corrections to C„ and M„ in (p^). 

The first major correction comes from the fact that we have taken into 
account in ( p3D that levels with rir + 1,1 = and rirJ = 2 are degenerate 
(in the lowest approximation when one neglects spin-dependent force) and 
hence used the sum of coefficients (0) for C„ in (^). However it is not true 
for the p-meson, where coefficient with / = yields | of C„. Therefore one 
should use instead of ip{z) in (P^D, ( P?]) the corrected expression 

^{z) tP{z) = iP{z) + ^—^ (100) 
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where ^—^ cancels the corresponding pole (with = 0) in (^) and replaces 
it with the correct coefficient 7 = 70 = I- 



Moreover, in the analysis it was realized that the value 70 = 3 (which 
does not take into account radiative correction to the p width) yields some 
10% discrepancy in the leptonic width of p meson. Now the destiny of the 
^ term in the asymptotic expansion (|9^) depends on the exact value of 7. 
Indeed, using (P^), one obtains the expansion 

=ln^+^(l-7-^) + 0(^) (101) 



? ^2 



•72- 

l2\ 



where z = 5-^ 

One can see that the requirement of the absence oi 0{1/Q^) terms implies 
the following condition of 7: 

2 1 , , 

7 = 3«:rad = - (102) 

where Krad = 1 ^i^, and An is the radiative correction due to gluon ex- 
changes in the leptonic width of p-meson. It is clear that condition is in the 
correct ballpark and one should take into account other possible corrections, 
to be discussed below. 

Till now we have considered the lowest approximation for for and 
Cn, Eqs.(^), (pTD , where is linear in n and C — n is n -independent. 



However, there are corrections to this quasi classical behaviour yielding 

= Anan + + ^ (103) 



a. ^ cr (1 + ;^) 



where C^^ is given in (|9T|). To see the effect of these corrections we ffist put 
b = and consider the sum (cf. Eq.(^)), which can be expressed through 
ip{z) again 

00 (1 + 3l) 

T. n. = 



n=nQ- 
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Thus we see that the correction ^ yields a power correction ^^""^^2 which 
modifies asymptotic expansion of ip{z) ( pH]) and can occasionally diminish 
some term to bring them closer in agreement with in the OPE (^). The same 
is true also for the mass correction b/ n in ( |103| ), which yields power correction 
(m^+q2)2 as a factor in expression like ( |104|) . However since coefficients in 
(p8| ) are fast diminishing, the overall cancellation in all terms seems to be 
unprobable and therefore the problem ii) remains unsolved. 



6 Analytic properties of ^^(s) 

In this section we shall discuss the analytic structure of the BPT in the 
limit Nc ^ oo m the whole complex plane of —s = and establish con- 
nection between our results and SPT. In what follows the basic role will be 
played by the function Il{Q^) and we shall use this example as a typical one, 
generalizing to other amplitudes at the end of this section. 

Our starting point is the BPT for n(Q^), Eg. (|8^) . Analytic properties of 
n*^°^(Q^) are given in Eq.(^), showing a sequence of purely nonperturbative 
poles 

s = -Q2 = (Mf)2,n = 0,1,2,... (105) 
Inclusion of to all orders makes the following changes with the un- 



perturbed spectrum (|105|) . The gluon exchanges play two different roles: a) 
instantaneous Coulomb interaction shifts NP poles b) gluons propagat- 
ing in the confining film form hybrid states. Correspondingly the NP poles 
are mixed and shifted by the hybrid states. Therefore one has the following 
hierarchy of the shifts of original NP meson masses 

M^°) M^°)(0(K)^)) ^ M^^ik) (106) 

The one-gluon hybrid state (apart from Coulomb shift, which we disregard 
for simplicity at this moment) with mass M^^^ is also shifted due to mixing 
with meson and higher hybrid states, the mixing being characterized by the 
power of a^, and therefore one has 

Mi'\0{{asf)) ^ Mi'\k) (107) 

In a similar way one has shifted states [k] for m-hybrid configurations 
(with m gluons). 
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As it is clear from the schematic example of section 4, the terms 0((?q) 
(the second term on the r.h.s.of (|57|)) refer to 1-gluon hybrid states and do not 
contain any analytic structure of asiQ'^), since they are not yet renormalized 
(and of course not yet RG improved, i.e.resummed). 

This renormalization process, leading to the nontrivial dependence of 
Q^b(Q^)) starts with the diagrams of the order 0{g^) in SPT (and of the 
order 0{V^V2) in the example of section 4 - see the third term on the r.h.s. 
of (57)). 

This term and its QCD analog has singularities of the double hybrid (the 
term Go{E — ui —002) in (57)), corresponding to the gluon-loop diagram and 
accompanying diagrams in OIqq) in SPT. 

Next one can do a RG improvement of this OIqq) result, which leads to 
the resummation presented in G^^^E), Eq. (68). In this approximation one 
can write 

47r 

as = = (108) 

6o[lnf + V^(g2)] ^ ^ 

where ipiQ"^) is a meromorphic function with poles at the double-hybrid po- 
sition, i.e. at 

s = Ml^^\k = 0) (109) 

In the particular case when ijj{Q'^) = ip (^~t^) and Mq = M^\k = 0) the 
form ( |108D coincides with the one, suggested in [|1^ . 
In the large limit one has 

and Ob has the same freezing form, which was extensively used before in 



12- 16 



47r 

= . 1 Q^+M,f (111) 
60 In ^2 " 

Now we can compare the form (|108D with our model solution (^6|). Indeed 



when is large, one can keep in the denominator of (^) the first two 

terms, and one has the form coinciding with ( |108| ). In the next approxima- 
tion one can keep in ( p.08| ) the term (^-^ + which yields next order 

approximation, when Mj^\0) is shifted to the higher order position M^\2) 
and so on. 
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Thus one can see that while the OIqq) term contains the double hy- 
brid poles, the RG improvement leads to the shifted poles, which do not 
correspond to the final physical poles, the latter are the result of infinite 
resurmmation to all orders. 

Let us now discuss the imaginary part of n(Q2). It is clear that U^^\Q^) 
contains the NP poles. The term OIqq) (equivalent of the second term on the 
r.h.s. of (|57|) ) contains double NP poles from {Go{E)Y and 1-gluon hybrid 
poles from Go(-E' — a;i). The O^Qq) term contains double NP poles, double 
1-gluon hybrid poles and 2-gluon hybrid poles - the latter are associated with 
ttfi in the lowest order. The same type of classification goes on for higher 
terms. The outcome for the expansion (|8^) is that in the term a'^U^'^^Q'^) 
part of singularities is associated with and another part with U^'^^Q'^). 

Thus we see that the requirements of i) logarithmic behaviour of ip{Q'^) 
at large and ii) meromorphic analytic properties in the whole plane, 
leads to the one possible representation of ip{Q^)] 

(112) 

where N is finite, so that logarithmic asymptotics of ip{z) is not modified. 

Let us now compare our Eq.(|3TD with the SPT, more explicitly with Gell- 
Mann-Low equation, written in the form adopted by Radyushkin 0. 




i = i„«! = i^ + !^,„^ + A + MLzMf!i + off^n, (113) 

A2 boas hi A-K 2hl Air \\AitJ J 

Here A is a parameter of integration which fixed the definition of A. This 
should be compared with our Eq.(^) which we can write in a similar form: 

, m2 /n2 + M2\ 47r 26i / 1 \ , , 



It is clear that our analysis with the expansion (|32D is equivalent to the 
analysis of (|113| ) for large and small a^. As we have found in (|39| ) the 
predicted values of 62*''^ , b^s^^ depend on the value of A and keeping the choice 
(^J-^ = ^ |^J_ _|_ a), one has from (||) for A > Aq = 1.145 

6^*''^ < 238.15 : b'i'"^ > 18266 (115) 
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To repair obtained values and get agreement with MS calculated values 
| 21| one can use the modified ip{z) function like that in (|112|) , namely 



1 \ , / 1 . \ A r 



\a.J Va, / r^^^ 1 + a^A + fc) 



fc=0 



It is important that the position of poles of ipik) are not changed, the only 
change is in the coefficients of those pole terms, which can change the ex- 
pansion (p7|), yielding agreement with MS coefficients of expansion (p^). 

We consider now the problem of analytic continuation of q;s(Q^) and 
q;b(Q^) from the Euclidean region > to the time-like region < 0. 

One can see that the freezing form ( |1 1 1| ) is well defined for > — Mg , 
and has a logarithmic branch point at = —Mq. 

However this is an artefact of representation ( p.l2| ), which is valid in the 

asymptotics when — ^3-^ ^ 1, otherwise one should use the original form 



dH). 

Choosing it in the form 

An 

one can see that ^^(Q^) is defined in all complex plane and has only 
isolated poles there at zeros of the denominator in ( p.l7|) . These poles are the 
result of partial resummation of two-gluon hybrid poles present in %() ( — ^2 ^ ) 
and lying at s = —Q^ = Mq -|- nm^, n = 0, 1, 2, ... We note that this resum- 
mation is the standard RG improvement of perturbative series, where large 
logarithms are summed up in the geometric series (or else a result of solution 
of RG equations). The crucial point is that most important physical thresh- 
olds are contained in 11*^°^ (Q^) and have nothing to do with singularities of 
asiQ^) (orof a,(g2)). 

In a similar way the singularities which appear from the asymptotic solu- 
tion of (|46|) for small as, namely term InaB(Q^) does not produce additional 
logarithmic cuts and is a result of asymptotic expansion of Eg. (^Gf) , where in 
exact solution only meromorphic function q;_b(s) is obtained. 

Thus as we discussed above the equation (HBI) or its generalized form 



An 

«B = -^—T. TTTT (118) 



bo [inf + ^(Q2) + ^(^) 
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where 'ip{Q'^) is given in ( |112| ) and (f is given in ( |116|) , produce ^^(Q^) 
which is meromorphic function of in the whole complex plane of Q^. 

Let us now look more carefully at the creation of imaginary part of 11 (Q^), 
which leads to the positive hadronic ratio R{s). We start with Yi'^'^\Q'^) (|85[) 
and (BBI): li^^^ in SPT and BPT coincide asymptotically, namely 



n(o)(^PT) = -^ln^ 



nW(i?PT) = -^In^-tiS (119) 
Defining s = —Q^, arg{MQ — s) = for s < Mq, one has for large s 

and for SPT one has on the r.h.s. instead of 0(s — Mq). 

Note however that the form (119) for 11*^°^ (BPT) is the asymptotics at 
large positive and in getting the result ( |12CI|) the analytic continuation of 
the logarithmic asymptotics was done unlawfully strictly speaking. 

To calculate in BPT more rigorously, one should use ( pSf ) in the region 
s > Mq, and one gets 



127r2 ^ \ m? I 127r2 - s + km? 



12vr ,=Q 



J2 5{Ml + krn^ - s) (121) 



One can introduce the average of /mn*^°^ = [s + ie] over some energy interval 
(this is the "duality interval" discussed e.g. in |^) comprising poles 

(/mn<»)(,„ + = £ ^-i i = ^ (122) 

One can see that (|122|) coincides with (|120D thus justifying the procedure of 



direct analytic continuation of the logarithms in ( |119| ), when the averaging 
over the duality interval is done. 

One can also see that the standard SPT expression ( |119| ) yields the correct 
answer for the same reason: the imaginary part of asymptotics (both of SPT 
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and BPT) coincides with the correct averaged imaginary part. We note that 
for SPT this check is impossible, since the correct procedure is not available. 

We now turn to the more delicate point of analytic continuation of a{s). 
In BPT, as was discussed above, as(s)is given by ( |118| ) or in "one-loop 



approximation" by (|117|) and is a meromorphic function, hence its analytic 
continuation is a direct and unique procedure. E.g. for ( p.l7| ) separating 
one pole in ijj (^ '^ ^ , = Mq + nm^, one has in its vicinity 

47r 

where Zq is to be found from the equation In ^ + ^p{zQ) = 0, and z = 



A2 

(If one uses exact model expression (^) instead of "one-loop approximation" 
Eq. ( |TT^ ) the r.h.s. of ([T23| ) and subsequent (Eqs. ([T2^ ) and ([^) should be 
multiplied by the constant factor (1 — ^ip'{A)), which effectively decreases 
ImaB)- Using relations from Appendix, one gets 

/.na.(.^M- + fe)^^;"'/'"^"-'-^°'"'' (124) 



y-^^ + SM + C(2) 
where 5q and Skiji) are defined in Appendix, 

2 

5^\n)=\n—~C + S,{n) 

Let us now compare this exact procedure with the suggestive one, using the 
asymptotic form of ( |117| ), Eq.( |lll| ). As before, we have 



asis + is) - asis - is) Att'^Q{s - M"^) 

Oojln I 1 + vr^l 

One can notice on the other hand, that averaging over some duality in- 
terval in ( |124|) yields 

47rM ^ 1 
bo A 6 ^{n + k) + ^ 

and this result as shown in Appendix asymptotically agrees with ( p.25| ), specif- 
ically at large s when the denominator in (|125|) grows. Now we turn to 
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another definition of discontinuity od wliicli is widely used in e~^e an- 
niliilation and r - decay. To make contact witli SPT calculations consider 
instead of U{Q^) the so called Adler function D{Q^) = ^^^^ with the 
expansion 



IT 



TT 



TT 



...} (127) 



where coefficient was found in pl|, and d2 in p8 



d2{MS) = 1.986 - O.llSn/, ds^MS) = 18.244 - 4.216n/ + 0.086nJ (128) 

The procedure suggested in to go from Di^Q"^) to R{s) is straightfor- 
ward, provided one knows discontinuity of n((5^) (see discussion in 0, 
and recent publication [p9l ) 

R{s) = — / D{a)- (129) 

where integration goes on both sides of the cut on real axis starting at s = 0. 

To proceed one inserts in ( |129| ) one - or two- loop expression for as{Q'^) 
and as the result one obtains an expansion for R{s) 



^(^) = Ee^ i + E4$ R (130) 




where ${} is the transformation defined by ( 129 ). In this procedure the InQ^ 
in the denominator of as(<5^) transforms into Ins with coefficients less than 
unity and asymptotically tending to one. 

Let us now consider what happens with our expression for asiQ"^) ( |117D , 
which is valid not only for large but in the whole plane. We assume 
that in BPT the function Db{Q^) is given by the same expansion ( |127D , 
where only should be replaced by asiQ'^) and we use ( |117|) for asiQ"^) to 



insert into ( |129| ). In this way one obtains, keeping the notation of 

ds' 




[^'{z^){z-z^)Ys'- 



;i3i) 
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As shown Appendix, iP'{zq) ~ (In MQ+^f]]_y_ Here the pole of the k-th 
order is at s' = Mq +nm'^ — 6Q{n)m'^ , and n satisfies condition MQ+nm? < s. 
For arbitrary k one obtains 

This is a property of the integral in ( |131| ) with the weight ^ that any function 
representable as a finite sum of poles yields zero since -^^ = 0. 

However in the asymptotic region of s = s(l + ^7), s ^ m^,when number 
of poles is large, one can use asymptotic form of ag, (|112|) , and the result is 

This result can be most easily obtained from (|129|), where the contour 



of integration Cr is modified to be a circle of radius R = s with center at 
cr = 0, namely 

R(s)^±l D(./A *((^)n = ^/ (134) 

The result ( |132D is obtained from ( |134| ) trivially introducing angular variable 
a = sexp{i(j)). For large s, s m^, (Mq)^ one can use in ( [134| ) asymptotic 
form of asicr) from (|111D . 
Thus one obtains for R{s): 



Ris) = N,Y.e 



1 



(135) 

Comparing with $|(^^^ | from Radyushkin one can see very close 



correspondence with our $ <{ ( — ~ 



7 Discussion and Conclusions 

We have developed in this paper the background perturbation theory (BPT) 
in the large limit, when all physical amplitudes and their perturbative 
expansions have only isolated singularities (poles). 
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Using this fact, we have proposed a nonperturbative solution to P{a) and 
ct{Q'^) containing all terms of loop expansion. 

By choosing a specific meromorphic function we have modelled both 
(3{aB) and asiQ^) satisfying necessary criteria of analyticity of in the 
Euclidean region of Q^, and analyticity of P{a) for positive values of as. 

We have constructed the time-dependent old-fashioned perturbation the- 
ory to investigate analytic properties of different terms of perturbation scries, 
and have shown that those reduce to isolated poles (of first and higher degree) 
shifted with respect to unperturbed bound state positions. 

This enables us to modify original ansatz for (3 (a) and asiQ^) to achieve 
a full correspondence with standard perturbative expansion of (3{a). 

Finally, we have studied the problem of analyticity of asiQ"^) and of ana- 
lytic continuation into region of time-like Q^, where asiQ"^) has singularities, 
and have compared this procedure with analytic continuation in SPT. 

We have found a striking similarity in the perturbative series for R{s) 
in case of SPT and BPT at large s. However our BPT series has correct 
physical singularities at large in contrast to SPT. 

More detailed calculations and comparison with experiment may decide 
on the applicability of both methods and are planned for the subsequent 
publication. 

Financial support of the RFFI grant 00-02-17836 and INTAS grants 00- 
110 and 00366 is gratefully acknowledged. 
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Appendix 



To find poles of q;b(s) 
the equation 



An 



bo 



m 



one must find the roots ZQ{n) of 



In — + i;{zo) = 



which can be represented as 

zoin) = -n + S, (5 < 1, n = 0,1,2. 
The following relations will be useful 

n 1 



k=0 



where Sk{n) = I]fc=i C = 0.577, and Eq.( |A.l|) takes the form 

1 m'^ 
---C + S,{n) + 5{1 + S2{n) + C(2)) = - In — . 



6 

To the lowest order one has 
S = So = 



A2 



+ 0(61). 



In a similar way one obtains iP'{zq), 

ij'{zo) = ^ + S2{n) + C{2) + 0{5). 



For n > l,S2{n) ^ C(2) 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



^(l + 5) = -C + 5{C(2)-5C(3) + ...} 

where C(2) = ^ , C(3) = 1.202. 

Hence one can write an expansion in powers of 6 

^{-n + 5) = -1 - C + 5i(n) + 5(1 + 82(71) + C(2)) + 0(5^) (A.5) 



(A.6) 



(A.7) 



(A.8) 
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In this way as in the vicinity of the pole is 

4:71 

bo[5o'\n) + S^in) + C(2)](M2 - S^m^ - s) 

where = Mq + nw? , 



aB{s ^ Ml + le) = , ^ , o , ...M.... F— ^ IV (A-9) 



Now for Si{n) the following asymptotics holds true 



S^{n) = C + \Yin + ^ + o(^]. (A.IO) 



Therefore ( |A.7| ) can be rewritten as 



^1 1^1 

" Inf + ln(n+l) = I^^Zpl " ^ ^"^'^^^ 

and the averaged imaginary part of in ( |126| ) has the form 

47r2 " 1 47r^ 

&oiV.^o(^ln^) +f 6o(ln^) 

where ~ s is in the middle of the averaging interval, which is asymptot- 
ically close to the expression in Eq. ( |125| ). 
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